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In this paper we show that a system of three fermions is exactly solvable for the 
case of a single- j in the presence of an angular momentum-J pairing interaction. On 
the basis of the solutions for this system, we obtain new sum rules for six-j symbols. 
When the Hamiltonian contains only an interaction between pairs of fermions coupled to 
spin J = J max = 2j — 1, the "non- integer" eigenvalues of three fermions with angular 
momentum I around the maximum appear as "non-integer" eigenvalues of four fermions 
if / is around (or larger than) J max - This pattern is also found in five and six fermion 
systems. A boson system with spin I exhibits a similar pattern. 
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1 Introduction 



To obtain a simple solution of the many-body Schrodinger equation is a long dream 
of physicists. There have been numerous efforts to obtain analytical solutions which 
do not require diagonalization of the secular equation by using computers. Among 
the many efforts in nuclear physics, the Elliott Model [1], the seniority scheme [2], 
the s and d interacting boson model [3] and a similar model by using schematic S 
and D pairs, the fermion dynamical symmetry model [4] are successful examples 
along this line. 

In Ref. [5] we showed that for a large array of states of four fermions in a single- j 
shell with only the J = J max pairing interaction, the eigenvalues are asymptotically 
integers labeled by numbers of J = J max pairs in the wavefunction, and that those 
corresponding wavefunctions of these low (I^ is the total angular momentum for 
a state of the four fermions) states are readily constructed in the nucleon pair basis. 
Besides the "integer" eigenvalues (as explained in Sec. 2), there are eigenvalues not 
close to integers when is around or larger than 2j, and little was known about 
these states. It would be desirable to discuss both the "integer" and "non-integer" 
eigenvalues on a more general footing. This is one of the goals in this paper. 

In Sec. II of this paper, we shall first study n = 3 systems (n is the number 
of fermions), which are readily solvable for any J pairing interaction only. The 
solutions of n = 3 provide an appropriate platform to explain the main idea of this 
paper. In Sec. Ill, we report relations between the eigenvalues of n — 3-, 4-, 5-, and 
6-particle systems with the J max pairing interaction. Using these relations one may 
obtain approximate values for both the eigenvalues and wave functions of low-lying 
states of these systems. We propose a hypothesis by which one readily obtains, to a 
high precision and without diagonalization of the Hamiltonian, the wavefunctions of 
states corresponding to some non-integer eigenvalues discussed in our earlier work 
[5]. A summary and discussion is given in Sec. IV. In the Appendix, we present a 
number of new sum rules for six-j symbols, some of which were derived recently by 
Talmi [6]. 

In this paper we use a convention that j (f) is a half integer, and that I (/') is 
an integer. They correspond to the angular momenta of the single-particle levels of 
fermions and spin carried by boson, respectively. J is used as the angular momentum 
coupled by two fermions in a single- j shell or two bosons with spin /; the maximum 
of J, J max , = 2j — I for fermions and 21 for two bosons. We use superscript (n) 
to specify the particle number n in the angular momentum 1^ and the eigenvalue 
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2 Three fermions in the presence of Hj only 

The pairing interaction which couples two fermions to an angular momentum J is 



Hj — Gj A M A J M , A M — —j= 



a) x a) 



j 



M=-J 

M = ~ (~ ^) M ~/K i^j X %1-M ' A J = J= [d,j X Clj] 



V2 



V2 



(1) 



where [ ] J M means coupled to angular momentum J and projection M. We take 
Gj = — 1 in this paper. 

We now consider the pair basis of three nucleons 



\ 3 *\jJ\lV,M) = -—(a]xA<% |0>, 



N 3%3J 



(2) 



where N^-j is the diagonal matrix element of the normalization matrix 



<,,= (0| (a, x 



(3) 



In general this basis is over complete and the normalization matrix may have zero 
eigenvalues for a given I^ 3 \ 



We first rewrite the matrix elements of Hj and Njj,.jj as follows [7] 
(f[ 3 K']I^,M\Hj\f[jK}I^,M) 



\f N jK';jK'\/ N jK;jK 

<0|( 



EH 



/(3)+J-L 



2/( 3 ) + 1 



X 



r(3) 1 L r rO^^A ^ 



|0>, 



7(3) 



A J ', (a} x A J t) 



/(3) 



|0>, 



(4) 



where L is a short hand notation of \/2L + 1. 

According to Eq. (12a) of Ref. [8] and Eqs. (3.10a) and (3.10b) [9], 



dj, (a'j x Qjj 



J 



A\ A s] 



rS TyS 5 tfi - 2fs 



r s t 
3 3 3 



(aj x djj 



(5) 



Using commutators in (5), we obtain 

L 



(0|[(s,x^') /(3) ,^tj =(0 | 

-2k>jj2(-y +K ' +L+J i" i3) i 



( \i<- 3) -jx x /(3) ~ 
J 

k' i& ) r k' j t 



J L t 



(-) l(3) - jI -^5 L , j (6 K , ;J + 2K'j 



K' J J 
J j 7( 3 ) 



J J J 



"3 ■ 



dj X 



(L)' 



(6) 



where a sum rule 

E(-) t+J+/(3) (2t + i 

t 

is used. The Hermitian conjugate of Eq. (6) yields 



K' J t 

3 3 3 



K' J t 

J J / (3) 



K' J J 
J j 1^ 



(7) 



A J , (a] x A*t) 



/(3) 



|o> = (-) 



2/(3) + j_ L . 



7(3) 

J 



# j j 
J j 7( 3 ) 



aj|0>. 
(8) 



Substituting Eqs. (6) and (8) into Eq. (4), we obtain 



NffL = + 2JJ 



J J J / (3) 
J' J J 



(f\jK']I^,M\Hj\f[jK]I^,M) 



\ LS jK';jK' 1S jK;jK 

Below we explain that there is at most one non-zero eigenvalue for each 7® of 
n = 3 with 77 = 77j. For a fixed J and for any 7®, we construct the \j 3 J : 1^)) and 
other states |j 3 TT : I^)) (K ^ J) which are orthogonal to |j 3 J : 1^}) as follows. 



\fJ ■ / (3) )) = \f\jJ]I {3) ) , 



N. 



1(3) 



\fK : 7( 3 ))) = |j 3 [^]7(3)) - ===== {K * J). (10) 

Using Eq. (9), we easily confirm that all matrix elements of the Hamiltonian in 
the basis (10), ((j 3 TT' : I^\Hj\fK : 7^)), are zero except ((\f\jJ]I {3) \Hj\j 3 [j J]I {3) )) 
= Njj^j. Thus all the eigenvalues of n = 3 for a given 7^ 3 ^ are zero for 77 = 77j 
except for the state with one pair with spin J, with an eigenvalue ^ (the 
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number in superscript specify the particle number n) given by —Njj.^j. This result 
was also proved recently by Talmi in terms of coefficients of fractional parentage [6] . 

Next we explain why the eigenvalues for the n = 3 cases are close to integers 
when H = Hj max . As shown above, the wave function of the lowest energy state 

for each /( 3 ) is given by ^(aj x aj^ max ' ) x afj |0). The eigenvalue Ej% J( .^ equals 

to —1 subtracted by a six-j symbol (refer to Eq. (9)), this six-j symbol is in fact 
very close to zero when is not close to I^l x = 3j — 3. The lowest eigenvalue for 
each 1^ (J® > j — 1) is thus very close to —1 unless ~ I^lx- All eigenvalues 
for /( 3 ) < j — 2 are zero. To show that the six-j symbols involved in Eq. (9) 
asymptotically vanishes, we list a few formulas of these six-j symbols: 

1 2j(2j-l)! 




I (4j - 1)! 

j(4j-3)(2j-i; 



(4j - 1)! 

j(4j 2 -4j-l)(2j-l)! 
(4j-l)! 



(11) 



which are less than 10 14 in magnitude for j = 31/2. Clearly, the approximate 
integer eigenvalues of n = 3 with H = i/j max comes from the fact that the six-j 

symbol < . . > are neglig ible unless J (3) ~ /® x - 

{ J J ^3 1 J 

The "non-integer" eigenvalues with ~ are also readily obtained: 



.tfW - 9 . 3 



4 3 L<W(j) 4 4(4j-3)' 

_ E (3) _ 27 _ _^5 21 

4 3 2x-2,J ma x(j) 16 32(4j-5) 32(4j-3) : 



£-(3) = _ , _ 

4 3 2x-3,j max (i) 16 32(4j-5) 32(4.7-3)' 



9 15 45 



-E {3) - _1 j 



81 105 15 1155 



64 512(4j- 7) 256(4j-5) 512(4j-3)' 
27 105 105 819 



.f] yo ' — — I 



(3) 

4 3 a l-5,J max (j) 32 256(4j-7) 128(4 j - 5) ' 256(4j-3)' 
279 ^ 315 ' 1785 



-E {3) = — - — + 

4 3 2x-6,J m ax(i) 256 4096(4j- 9) 4096(4j - 7) 

9135 17325 



4096(4j-5) 4096(4j-3)' 
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(3) 243 945 315 



i£L-7,Jm**(j) 256 4096(4j- 9) 4096(4j - 7) 
17325 21879 



4096(4j- 5) 4096(4j-3)' 

(3) 1053 3465 12915 

4Xx-8,J m ax(j) ~~ 1024 + 131072(4j - 11) ~ 32768 (4j - 9) 

58905 225225 855855 

+ 



65536(4,7- 7) 32768 (4j - 5) 131072(4; - 3) ' 



_ (3) 63 3465 3465 



/^x-9,J m ax(j) 64 32768(4j - 11) 8192(4j - 9) 

45045 83655 255255 

+ 16384(4j - 7) ~ 8192(4j - 5) + 32768(4j - 3)' ^ ' 

etc. We see that these above eigenvalues E^h , ,., stagger and saturate at —1 as 

J- v ; ) <J max \J ) 

1^ becomes smaller and smaller (but > j — 1). In the large j limit, the non-zero 
eigenvalue for each 1^ takes the first term; for a very small j value, e.g., j = 9/2, the 
eigenvalue %. Vm(9/2) (i.e., I® = 7/2, £$ )Jm(Bc(9/2) equals to -gf) is already 
very close to —1 (within a precision of 10~ 2 ). This explains why we frequently 
obtain asymptotic —1 eigenvalues for H = -£/j max and n — 3. For a state which has 

(3) 

_E r(3) . , ~ —1, the corresponding wavefunction can be understood as a single-?' 
"spectator" coupled to one pair with spin J max - There is no such a spectator for 
~ i£2x states, although their wave functions can be written as j 2 (J ma , x )j : J®). 
We note that three bosons with spin / exhibit a similar pattern: there is up to 
one non-integer eigenvalue for each I in the presence of boson Hamiltonian Hj, 

i$> >J(I) = -l-2(2J+l){ J \ ^ } , (13) 



lven by a pair witri spin J , i.e., \ b\ x ty'j 

coupled with a single boson operator b\. 



i.e., this non-zero (zero) eigenstate are given by a pair with spin J, i.e., (bj x tfj 1 



3 Relations between states of n = 3 and 4 for H = 



#7 



max 



In this section, we first discuss the cases with n — 4. In Ref. [5], it was found that 
the eigenvalues of n = 4 are asymptotically 0, —1 or —2 for small I^. These states 
are constructed by coupling one or two pairs with spin J = J max . However, some 
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"non-integer" eigenvalues appear as 1^ is larger than 2j — 9. These values are very 
stable for 2j — 8 < 4j — 12, and the origin for these states was unknown. 

Let us compare the eigenvalues of a system with n = 3 and n = 4 fermions with 
H = Hj max for j = 31/2. The distribution of all non-zero eigenvalues for n = 3 and 
4 is plotted in Fig. l(a)-(c) 1 , where (a), (b) and (c) correspond to the range of \E\ 
from to 0.8, 0.8 to 1.5, 1.5 to 3.8, respectively. 

From Fig. 1, we see that these eigenvalues are clustered at a few values but with 
exceptions. The "clustered" values are very close to the eigenvalues of n = 3 2 . This 
indicates that the eigenstates of n = 4 are closely related to those of n = 3. 

For j = 31/2 and n = 4 the total number of states is 790. The number of states 
with non-zero eigenvalues is 380. Within a precision 10~ 2 , the eigenvalues of these 
308 states are located at the eigenvalues of n — 3, and 21 states have eigenvalues 
closely at —2. We note that almost all the "non-integer" eigenvalues of n = 4 can 
be rather accurately given by one of three-particle clusters with 1^ ~ I^l x coupled 
to a single-j particle. In this example only four states with 1^ = 48, two states 
with /( 4 ) = 46, and two states with 1^ = 44 cannot be understood by either one 
of three-particle clusters with 1^ ~ J^x coupled to a single-j particle or two pairs 
with one or two spins being J max - 

For example, the peak for n = 4 in Fig. 1(c) near 2.25 is very close to the 
energy of {E^ | of n — 3. For j = 31/2, the maximum angular momentum 

-^max > J max {j ) 

I® of three fermions is f . The E { % =-ff|=-2. 26271186440677966. The 

2 4 3 ax,Jmax(i) 118 

minimum 1^ obtained by coupling a three-body cluster with = /® x t° a single- 
j particle is given by 3j — 3 — j = 2j — 3 (the triangle relation for vector couplings) 
and here 28. We find that the lowest eigenvalue of 1^ = 28 for n = 4 obtained 
from a shell model diagonalization is -2.26271186440689. The E% (close 

to -2.26) with J( 4 ) between 28 to 56 are listed in Table I. Two observations can be 
made: (1) the lowest state of each 1^ are well separated from the second lowest 
one, and (2), there is no eigenvalue which is smaller than -2.0 when is lower 
than 2j — 3 for n = 4 3 . 

We also see that the overlap of the wave function obtained by the exact shell 

1 The inset in Fig. 1(b) is rc-scalcd in order to see more clearly the exceptions of energies which 
are not close to those of n = 3. 

2 There is one peak at 2.0 which was explained using two pairs with spin J = 2j — 1 in Ref. [5]. 

3 We also note that the above nearly equality is asymptotic for a rather large j, not exact. For 
examples, when j is very small, for j = 7/2, the energy of ^ is — 1^=2. 31818182 while the 

lowest energies of 1^ = 2j — 3=4 for four fermions obtained by diagonalization is — | = 2.66666667; 
for j = 9/2, the energy of E^l } is — 2| while the lowest energies of = 2j — 3=6 for four 

fermions obtained by diagonalization is —2.34965034965038. 
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model with the state constructed as the /® x state coupled to a single-j particle, 
[Imlx x j] /(4) ; is ver y close to I. This may be argued as follows. The eigenvalue for 
the lowest state of n = 4 for 1^ > 28 is very close to the matrix element of the 
Hamiltonian for the state [I^l x x j] /4> . Suppose that the lowest spin 1^ state is 
not degenerate. Then a certain state which produces the same energy as the lowest 
spin /( 4 ) state will have the same wavefunction. In Table I most of the energies 
obtained by diagonalizing H Jingjc for n = 4 are close to the matrix element for the 
pure configuration of the [I^l x x j] /(4> state (also close to E^s) ). Overlaps of states 
having other "non-integer" eigenvalues near —2.25 for n = 4 with those given by the 
[/max x j] /(4> are close to 1, except three cases (two of them can be approximated by 
other three-particle clusters with 1^ ~ J r ^ x (but 1^ ^ J^) coupled to a single 
particle). 

We have calculated all overlaps between states of n = 4 which have energies 
close to the peaks and those of simple wavefunctions obtained by coupling a single 
particle to a non-zero energy cluster with ~ I^l x of three fermions. These show 
a similar situation as Table I. Therefore, we conclude that those stable "non-integer" 
eigenvalues of n = 4 with H = i^j max in Fig. 1 are given to a high precision by a 
three-particle cluster (nonzero energy) coupled to a single-j particle. 

One may ask which picture is more relevant to the states of n = 4 with eigen- 
values close to integers; one in which a three-particle cluster (nonzero energy) is 
coupled to a single-j particle, [1^ x j] /(4) with J® ~ I^Ixi as proposed in this pa- 
per, or one in which four particles are coupled pairwise with one or two spins being 
Jmax, as proposed in Ref. [5]? 

First we note that for > 53 only a single state is possible and these two 
pictures are therefore equivalent and exact; for = (or 3) the number of states 
is the largest integer not exceeding (2j + 3)/6 (or (2j — 3)/6) which is larger than 
1 in most cases [10] but there is only one state which gives a non-zero eigenvalue 
for the J max pairing interaction [5]. Also in this case the two pictures are therefore 
equivalent and exact. 

For states with 1^ < J max and energy close to —2.0, it was proven in Ref. [5] that 
a description by using two pairs with two spins being J max is very good. For these 
states one may ask whether a description by using a single-j particle coupled to one of 
three-particle clusters with 1^ ~ J r ^ x is also relevant. To see whether or not this is 
true, we calculate the overlaps of the states of four fermions with E^ 4) J ^ ~ —2.0 

and < J max which were obtained by the shell model diagonalization with all 
possible three-particle clusters coupled with a single-j particle. These overlaps are 
between around 0.6-0.8. Thus the picture using a single-j particle coupled to the 
three-particle cluster for states with 1^ < J max and energy close to —2.0 is not 
appropriate. 
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Then, how does the picture of two pairs with spin J = J max work as increases? 
We calculate the overlaps of states with energy around —2. We see that for 1^ = 42 
( = Imlx ~ 14) the overlap is still 0.9962, showing that the pair picture is still very 
good. 

The next question is related to the states of four fermions with energies near 
— 1. There are about 100 states for n = 4 and j = 31/2. Both pictures can give 
eigenvalues at —1. The number of states with E ~ — 1 is not unique. As was shown 
in Ref. [5], for states with small the number of states with E^l } J ~ — 1 

states is the largest integer not exceeding 1^/2, which is larger than 1 except for 
= 0, 2 and 3. Because these eigenvalues are very close to but not exactly — 1, 
the mixing of these configurations can be large. Coupling two pairs, one with spin 
Jmax and the other spin J' ^ J max picture gives a very good classification of states, 
but not the exact wavef unctions. On the other hand, the picture using a three- 
particle cluster of nonzero energy coupled with a single- j particle provides us with 
better wavefunction than the pair picture, but it does not provide us the number of 
states with E^ 4) } ... ~ — 1. These two pictures are therefore complementary in 

describing the states for n = 4 with e[% t , ., ~ —1. 

4 States of five particles and those of six particles 
with H = ff Jm „ 

In this section, we proceed to more particle systems. Although we did not find 
simple descriptions for them, we are able to find some relations between states of 
different n systems with an attractive J max pairing interaction. 

The picture using clusters of N {M < n) particles coupled to (n — A/") single-j 
is also found in states of systems with n > 4. We study in this section the j = 19/2 
shell for both n = 5 and 6. The cases with larger j-shells yield a similar picture 
with higher accuracy. 

Asymptotic integers appear in the eigenvalues e[% t , when /( 5 ) is not very 
large. They are either zero, or very close to —1 and —2. The number of states for 
1^ = | is three, among which there is one with zero eigenvalue, one with eigenvalues 
~ —1 (within a precision of 0.01) and one with eigenvalue ~ —2 (within a precision 
of 0.01). The number of states for = | is seven, among which there are two 
with zero eigenvalues, three with eigenvalues ~ —1 (within a precision of 0.01) and 
two with eigenvalues ~ —2 (within a precision of 0.01). A similar situation holds 
for larger 1^ states except that eigenvalues ~ (E$) j ^ ~ 1) (1^ ~ -^mlx) or 
£^3) j ^ appear. Corresponding to each three-body cluster with ~ I^l x , the 
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minimum of I which gives "non-integer" eigenvalues EfX . , ~ (E^h , ,., — 1) 
and ^3) iJmax0) ~ ^S) >Jmax0 -) is / {3) - (2j - 1) and - (2j - 3), respectively. For 
examples, E^L T , ., ~ £^ 3 i — 1 = ||f| = —1.89707 appears in states with 

"» * e-S-W-1) =i"7 = I: fig^o, ~ £|, _ 5 , Jm „ 0) = if = -0-89707 
appears in states with 1^ > /^x — 5 — (2j — 3) = j — 5 = §, etc. The non-zero 

(3) 

eigenvalues for n = 5 are equal to or concentrated around 0, — 1, —2, ~ -5/(3) j ,^ 
and ~ (-^j(3) j _ 1) with ~ 4^. The above regularities survive unless 

max ' 

Now let us look at the case with n = 6 in the same shell j = 19/2. Below we 
consider the case with = as an example because other low states behave 
similarly. There are ten states with = 0. Among them there are two with zero 
eigenvalues. Non-zero eigenvalues are more complicated than systems with smaller 
n, because n = 6 can be divided into more sets of clusters. We first divide n = 6 
into two clusters with n = 3 and ~ /£2x for each three-body cluster. Then 
we obtain eigenstates with eigenvalues: -4.54286, -3.31055, -1.23269, -2.42027, 

(3) 

— 1.75573, and —2.10284, which are approximately equal to twice those of E V I{ ' 3) 3 , 

(/ (3) - 41): - -f ,-f respectively. We canTlso 

divide 6 into three two-body pairs. Here we take 1^ = J max for these three pairs 
which lead to an eigenvalue very close to —3 (—3.01537). Besides these eigenvalues, 
there are one eigenvalues which are close to —1.0. 

We did not succeed in setting up a simple scenario of the distribution for all 
eigenvalues of systems with n = 5 or 6 and H = H Jin ^. This is partly because the 
number of states for each is not analytically known. The number of combinations 
for different clusters is also much larger than the cases with n = 3 and n — 4. 

Based on these relations we suggest that the low-lying states of each 1^ of 
fermions (bosons) in a single-j shell (with spin I) interacting by an attractive H = 
ifj max favor a cluster structure, where each cluster has a maximum (or close to 
maximum) angular momentum. The coupling between the constituent clusters (in- 
cluding pairs and spectators) are very weak and negligible, therefore we can obtain 
both their approximate wave functions and eigenvalues, which are simple summation 
of those of the clusters. 

In Fig. 2, we showed the distribution of all non-zero eigenvalues for systems with 
n ranging from 2 to 6. It is easy to notice that the eigenvalues are concentrated 
around some values for n = 2 to 5. This pattern becomes less striking for n — 6. 
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5 Discussion and summary 



In this paper, we first show that a system of three fermions in a single-j shell in the 
presence of H = Hj is solvable. We prove that there is at most one state with a 
non-zero eigenvalue for each . We can analytically construct both the eigenvalues 
and corresponding wave functions. A similar remark applies to three bosons with 
spin / in the presence of Hj. On the basis of the above results for n = 3 a series of 
new sum rules of six-j symbols can be found. 

We show that the eigenvalues of three fermions in a single-j shell with H = Hj max 
are very close to or —1 unless J® ~ I^l x = 3j — 3. This kind of situation is very 
similar to the case of n — 4, as studied in Ref. [5]. 

We also find that the "non-integer" eigenvalues of ~ for n = 3 appear 
as "non-integer" eigenvalues for n = 4 when is around or larger than J max . The 
overlaps between the wavefunction of these "non-integer" eigenvalues of n = 4 and 
that of ~ /£2x state coupled to a single-j particle is very close to 1. This finding 
allows us to construct approximately the states of n = 4 by using results of n = 3 
as we have shown. We confirmed that this is also true for five and six fermions in a 
single-j shell in the presence of J max pairing interaction. Bosons with spin / exhibit 
a similar pattern. Similar regularity was found for n = 5 and 6, although we did 
not succeed in setting up a simple rule for all states. 

The relations between E^L T c-n , ,.,e[% , indicate the follow- 

ing pattern: the attractive J max pairing interaction favors clusters (including pairs 
and spectators), where the angular momentum of each cluster is close to the max- 
imum. One thus explains the "integer" eigenvalues and "non-integer" eigenvalues 
proposed in Ref. [5] by using a picture of the clusters for fermions in a single-j shell 
or bosons with spin /. 

As is well known, the existence of degeneracy indicates that the Hamiltonian 
has a certain symmetry. The degeneracy for the J max pairing interaction, however, 
is not exact. It would be interesting to explore the broken symmetry hidden in 
the J max pairing interaction discussed in this paper. It would be also interesting to 
discuss the modification of the J max pairing interaction in order to recover the exact 
degeneracy. 

Acknowledgement: We wish to extend our special thanks to Prof. I. Talmi 
for his valuable comments concerning the sum rules of six-j symbols. We also thank 
Dr. N. Yoshinaga for discussions in the early stage of this work, and Drs. I. Talmi 
and O. Scholten for their reading of this manuscript. 
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Table I The lowest eigenvalues of the 1^ states in a single-j (j = 31/2) shell 
with /( 4 ) between 28 to I^l x = 56. When 1^ is smaller than 48 there is no eigenvalue 
lower than —2. The eigenvalue of the 7®x state with three fermions in the same 
single-j shell is -fff =-2.26271 186440677966. The column "(SM)" is obtained by 
a shell model diagonalization, and the column is matrix element of Hj max for 
the state constructed by three-fermion with J® = 7^x coupled to a spectator. The 
column "error" presents the difference between Ej and Si (two effective digits). The 
column "overlap" is the overlap between the lowest eigenstates of n = 4 and the 
states obtained by coupling single fermion a] to the I^l x state. Italic font is used for 
three cases for which the overlap is not close to 1. We note that the case of 1^ = 50 
(52) can be approximated rather accurately (10~ 3 ) as a three-particle cluster with 
= i^2x ~ 8 (/£2x — 6) coupled to a single-j spectator. 
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Ei (SM) 



Ej (coupled) 



"error" 



1.1 x 10~ 13 
4 x 10~ 14 

1 x 10" 14 
0.9 x 10" 13 
l.lxl0~ 13 

2 x 1(T 14 
1.9 x 10~ n 

8 x lir 14 

1.1 x 10" 9 
6.6 x 10~ 13 

4.3 x 1(T 8 
2.8 x lO" 11 

1.2 x iir 6 

9.6 x 10- 10 

2.4 x 1(T 5 

2.7 x 1(T 8 

4.0 x 10~ 4 
6.2 x 10~ 7 

6.1 x 10~ 3 

1.2 x 10" 5 
7.7 x 10~ 2 
1.1 x 10" 4 
2.6 x lO" 1 



3.0 x 10- 1 






overlap 



-2.26271186440689 
-2.26271186440682 
-2.26271186440678 
-2.26271186440669 
-2.26271186440692 
-2.26271186440700 
-2.26271186442899 
-2.26271186442233 
-2.26271186573172 
-2.26271186512903 
-2.26271191546116 
-2.26271188689249 
-2.26271325181426 
-2.26271236805292 
-2.26274016611845 
-2.26272031287460 
-2.26317530567842 
-2.26282037299297 
-2.26963309159052 
-2.26378385186917 
-2.34719850307215 
-2.27068252318197 
-2.57872583562800 
-2.30488200470359 
-2.89017281282010 
-2.41926851025870 
-3.24511394047522 
-3.66369313113292 



-2.262711864406782 
-2.262711864406777 
-2.262711864406780 
-2.262711864406782 
-2.262711864406805 
-2.262711864406981 
-2.262711864409884 
-2.262711864422405 
-2.262711864593695 
-2.262711865128374 
-2.262711871692375 
-2.262711886864667 
-2.262712064607266 
-2.262712367094411 
-2.262715960012236 
-2.262720286322401 
-2.262776481261782 
-2.262819747102017 
-2.263514816015588 
-2.263772302947436 
-2.270625142453812 
-2.270571840272616 
-2.323429204525185 
-2.304882004703592 
-2.592166600952603 
-2.419268510258698 
-3.245113940475225 
-3.663693131132918 



1.000000000000000 

1.000000000000000 

1.000000000000000 

0.999999999999999 

0.999999999999659 

1.000000000000001 

0.999999999963606 

0.999999999999996 

0.999999997833653 

0.999999999999758 

0.999999916591887 

0.999999999988818 

0.999997726566892 

0.999999999598199 

0.999952666087478 

0.999999987392690 

0.999151747579904 

0.999999632523561 

0.982828211942919 

0.999992036003522 

0.780582505446094 

0.999929221753443 

0. 706859839896674 

1.000000000000000 

0.873170713095796 

1.000000000000000 

1.000000000000000 

1.000000000000000 
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Appendix A New sum rules of six-j symbols 



The solution of Hj for n = 3 gives new sum rules. The procedure to obtain these 
sum rules is straightforward. As is well known, the summation of all eigenvalues 



with a fixed I is equal to 



n(n— 1) 



times the number of I states, where n is the particle 



number. For n = 3, the number of states can be empirically expressed in a compact 
formula [10]. 

In Ref. [5] we applied this idea and obtained that 



E (2-7+1) 

even J 



3 3 J 
3 3 J 



2j + 3 



6 



2j 



\ if 2j = 3k, 
if 2j" = 3A;+1, (14) 
-\ if 2j = 3fc + 2, 



where j is a half integer, and [x] means to take the largest integer not exceeding x. 
we derive a similar sum rule using the 1 = states of four bosons with spin I: 



E (2-7 + 1) 

even J 



j I I J ) 


3 


"f 




\l I J 




3 





f 1 ifZ = 3ib, 
| if Z = 3A; + 1, 
if I = 3k + 2, 

Below we give other sum rules of six-j symbols. For a half integer j, 

'21+3' 



(15) 



E 2(2J + 1) 

J=even 





I 




f3 




3 




I 3 



6 

3J-3-J 
6 



+ 35} 



3j+l-J 
2 



if 7 < j; 
if 7 > j . 



(16) 



where 



if (3j -3-1) mod 6 = 1 

1 otherwise . 



For integer /, we obtain similar sum rules given as follows. For I < I (I is an 
integer) , 



E 2(2J+1) 

J=even 



where 



I I J 
I I J 



I 

3 

3/-/ 



+ 35\ 



31-1+2 
2 



if (3/ - I) mod 6 = 1 

1 otherwise . 



if I < I; 
if / > I . 



(17) 



It is noted that the sum rules (14) and (15) are special cases of the sum rules 
(16) and (17). 
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Starting from Eq.(10.14) of Ref. [11] for J' = J", J x = J 3 = j, J 2 = J 4 = /, we 
multiply (2J + 1) and sum over J'. Using Eq. (10.13), we obtain 



£(2J + 1) 
j 



l V J 



>=0; 
= 1 . 



(18) 



Similarly, we obtain 



£(2J + 1) 
j 



I j J 
I j J 



and 



£2(2,7+1) 
j 



J I J 

j j J 



-1 if/<j, 
otherwise . 



if / < j , 

(-y-j-i if i> 3 . 



(19) 



(20) 



Using Eqs. (20) and (16), we obtain 



£ 2(2 J + l) 

J=odd 

and for I > j, 



{) j j} 



/+ 2- 3 



2j + 3 



' -1 if 2j = 3k , 

if 2j = 3k + 1 , (21) 

1 if 2j = 3A; + 2 . 



£ 2(2J+1){^ ^ J \ = (-y-i-l-3 



3J-3-I 



-35] + 



3j + 1 - I 



(22) 



Similar to Eq. (20), we obtain 
£2(2 J+l) 



I I J 
I I J 



2(-y+ l if I < I , 
1 + if I > I . 



(23) 



Using Eqs. (23) and (17), we obtain that 

:-y+ i + i-3 



£ 2(2J + 1) 

,/=odd 



I I J 
I I J 



(-)'* - 3 



3«-/ 



- 1 

35\ 



31-1 



if / < I ; 
+ 2 if / > / . 

(24) 
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Caption: 



Fig. 1 Detailed distribution of all non-zero eigenvalues for n — 4. The inset 
in Fig. 1(b) is rescaled to distinguish a few exceptional cases which energies are 
not close to those of n — 3. (a), (b) and (c) corresponds to different range of 
\E ( T t , ,_.J. We see that the eigenvalues for n = 4 are "clustered" at those of n = 3 
with few exceptions. 



Fig. 2 Distribution of non-zero eigenvalues \E^L , , .J for systems with n 
ranging from 2 to 6 and j = 19/2. One sees that the non-zero eigenvalues are highly 
concentrated. The concentration of eigenvalues for the case of n = 6 is less striking. 
The distribution is plotted using the number of counts for each J ^| (with 

the step length being 0.01) divided by the total number of non-zero eigenvalues. For 
H = Hj max with j = 19/2, the number of non-zero eigenvalues is 1, 17, 122, 472, 
1224 (in comparison with the number of the shell model space: 10, 45, 177, 521, 
1242) for n — 2, 3, • • •, 6, respectively. 



17 



Figure 1 Nov 20th/2003 



20 - 



(a) 



0.2 0.4 



ditribution of eigenvalues of n=4 



0.6 



j eigenvalues of n=3 



II 



CD 

_^ 

> 

CD 
D) 

'o 

CD 
CD 



10 - 



c 

O 

O 



0.8 

- 20 

- 10 




Eigenvalues of n=4 



